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ABSTRACT 

We have considered precession in accretion disks in which a second moment of 
inertia relative to an axis perpendicular to the axis of rotation may be very impor- 
tant. This formalism, that takes into account the precession contribution to the 
angular momentum, is based on the existence of a parameter p which determines 
three characteristic densities resulting from the averaging process and imposes 
constraints on the actual disk density. It is shown that the precession velocity 
will lie in a three branch solution, and depends on how large is the disk actual 
density as compared to the characteristic densities. Besides the large spread on 
the solution for the precession velocity, depending on the density strength, it 
may be prograde and retrograde. It is shown that the keplerian thin disk, with 
very large density values compared to characteristic ones , only precesses very 
far away from the primary object, which implies very large precession periods. 
For other models, the disk will thicken, with large deviations from the keplerian 
approximation. Constraints on the density only will be effective for very large 
values of the ratio -M-, respectively, accretion rate and mass of the primary. Un- 
der this condition, the structure of the precessing region is found. Lower bounds 
on the precession period are found for not so large values of this ratio. Deviations 
from the mean precessional motion are considered. It is shown that these devia- 
tions result in periodic motions as long as the time scales associated to them are 
comparable to the remaining time scales. Otherwise, they result in misalignment 
motions, forcing the plane of the disk to become normal to the orbital plane of 
the secondary. 



1. Introduction 



Precessional activity has long been invoked to explain time variations occurring in the 
spectra of galactic X-ray binaries like LMC X-4, Her X-l, SS 443, Cygnus X-l, etc. Obser- 
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vational evidences give support to the idea that these systems, and quite a lot of suspected 
others, consist of a disk-like system plus a third body. This is very suggestive of precession 
as being the effect due to the perturbing torque of a distant companion star in the the disk 
or, even, as the result in the disk of the secondary star precessing in the tidal field of the 
central compact objec t . As a matter of fact, the scenario is not unique, and according to 
Priedhorsky and Holtl (119871 ). there are, at least , five possibilities. This is so because there 
is uncertainty concerning the reason why the disk precesses, where and the extent on it that 
precesses. The possibilities explicitly dependent upon tidal torques under which the disk 
will precess can be summarized by the following two scenarios: 

a- the disk has a permanently tilted edge that freely precesses in the tidal field of the 
companion star. This kind of scenario has been first exploited by iKata (119731 ). tryi ng to 
explai n the 35 day period on the light curve of Her X-l. This same model was used by iKatz 
(119801 ) to explain observational features in SS 443; 



b- the companion star precesses in the tidal field of the central compa ct star, s l avishl y 
followed by the the outer tilted disk. This kind of scheme was proposed by [Roberta (119741 ). 
also to explain the 35 day period observed in Her X-l. 

However, precession is not a exclusiveness of galactic X -ray binaries sources, its presence 
being suspected in a lot of other astrophysical systems. Jet stru ctures obser ved in AGN, and 
also in protostars seem to be associated with a precessing disk ( iLivid (ll999l )). In the context 
of disks aro und supermassive objects, the reason why the disk precesses is still more unclear. 
Katzl (119971 ) suggested that observational data in O J 287 could be explained by a tidal torque 
due to the presence of another massiv e co mpanion star, i n a way similar to that employed 
to Her X-l and SS 443. iPringld (119961 ) and iPringld (119971 ) su ggests precession as a result of 
a warping ins tability cau s ed by the central sour ce irradiation. iRomero et all (120001 ) argue, in 
analogy with IKatz et all (119821 ) and iKatzl (119971 ).. aiming to explain the pattern of ejection 
observed in the Quasar 3C 273, that precession is due to a massiv e secondary black hole 
tidal t orque on the disk, a model, subsequently, applied to 3C 345 by lCaproni and Abraham 
(l2002h . 



It should be argued, however, that precession is a kind of solid-like response of the disk. 
A fluid-like response would be a more appropriate one for a fluid system like the disk. Be- 
sides, in the case of fluid-like response, the tidal torque acts as a perturbation, which main 
consequence is the generation of different sort of waves. Quite a variety of physical issues, 
suitable for a wave approach, such as several tidally driven instabilities, disk tidal deforma- 



tion (warpin g), horizontal 



and so on fdLubow 



1991 



Vishniac and Diamond 



l y and vertically driven resonances, angular momentum transport 



19921 ; iLubow and Pringld 11993 lOgilvie al 12002 ; lOgilvie bl 12002 



1992])) have been tackled quite successfully. In addition, it has been 
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shown that, under this approach, for disks in systems with extreme mass ratio, preces- 



(Lubow 


( 


1992 





that the axisymmetric part of the tidal potential is related to the disk solid-like response, 
while the disk fluid-like response is related to the non-symmetric part of the potential. In 
these works, the the disk solid body-like response comes as an integrabil;ity condition. 



Katz et all (119821 ) proposal to compute frequencies treats the disk as a ring, with no 
allowance given for the width in the plane of the ring, nor to the thickness perpendicular 
to that plane. By employing a suitable averaging procedure, he obtains for the for the 
precession angular velocity 



n 3 ^ 2 a 
Up = -- — cosd, 

4 "J 



where uo s is the keplerian angular velocity of the secondary, is the ring angular ve- 
locity, and 5 is the inclination angle. This formalism, besides treating the disk as a rigid 
body, not taking into account its structure, assumes the disk angular velocity parallel to the 
disk angular momentum . Though keeping fl p parallel to L, ass uming thin disks obeying 



polytr o pic equation of s tate, t he structure has been cons i dered by iPapaloizou and Terquem 



( 119951 ). lLarwood et al (|l996h . lLarwood and Papaloizoul (119971 ) and lLarwood! (119971 ). De 



composing the secondary tidal potential into odd and even z parts, with subsequent Fourier 
decomposition, and arguing that the odd axisymmetric part is responsible for rigid preces- 
sion, these authors, assuming - constant, obtain 



3 
4 



7- 2n 



n 



— cosd, 



as an integrability condition. In this expression, n is the polytropic index and fid is the 
disk angul a r velo city evaluated at the outer edge of the disk, R . For n — 2, they recover 
Katz et all (119821 ) result. According to them, the validity of this result holds whenever the 
disk doesn't thicken, the perturbative tidal potential is weak, and, above all, the sound 
crossing time is shorter than the precession period. It should be argued, however, that 
the disk angular momentum and the precession velocity are not parallel. Besides, heat 
generation, cooling and angular momentum transport do occur in accretion disks, a fact 
that makes not realistic the structure obtained from a polytropic equation of state. This 
imposes severe restrictions on the model of disk we want to consider.lt is reasonable to 
expect the polytropic index varying from n = |, for a (monoatomic) gas pressure dominated 
disk, to n = 3 for a radiation pressure dominated one. Disks with negligible temperature 
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gradients would require a polytropic index much greater than 3: 7 = 1 is appropriate for 
an isothermal atmosphere ( n — ► 00). In that situation, the mass content in the disk, the 
moment of inertia and the precession velocity are determined by values close to the inner 
radius rather than close to the outer radius, the disk structure being very susceptible to 
relativistic corrections and gravitational radiation losses may become important. Besides, if 
one is interested in a region in parameter space M — M where the disk thickens, the dynamics 
can not be dissociated from the structure, as is the case for the keplerian thin disk. In other 
words, the azimuthal velocity depends on the disk height scale . In addition, for a given tidal 
potential, its relative strength increases as the disk thickens. For actual radiation pressure 
dominated, the height scale is constant, the polytropic index is negative, implying negative 
specific heat. Besides this, in this case, one can not satisfy null boundary condition at the 
surface of the disk Finally, the precession period diminishes relative to the sound crossing 
time as the disk gets thick. A final possible question concerning the applicability of of these 
results to disks is related to the precession velocity dependence on the outer radius of the 
disk, which yields proportionality with the mass of the disk. Despite implying larger inertia 
moments, the larger the mass the faster the disk will precess. 



Romero et all (120001 ) work consists essentially of the application of iKatz et all (119821 ) and 



Katzl Jl997h to a binary system composed by a disk around a primary supermassive black 



hole, plus a supermassive secondary orbiting the primary in a keplerian way, not coplanar 
to the disk. It should be argued, however, that the formulation followed by these authors, 
besides suffering from the drawbacks, we just mentioned, of any precessing disk model, relies 
heavily on the above cited authors results for the precession angular velocity, which was 
obtained without considering the precession contribution to the total angular momentum. 
Besides, this model is supposed to apply to disks eventually dominated by radiation pressure. 
In addition, account of the precession contribution to the angular momentum introduces 
some non-linearity in the problem, which, in turn, introduces conditions on the averaging 
procedure, imposing constraints on the density, on the disk unperturbed angular velocity, on 
the viscosity parameter as well as on the height scale of the disk. In that case, there will be 
a contribution from the precession velocity to a motion in a plane perpendicular to the plane 
containing the normal to the disk and the normal to the plane of the orbit of the secondary. 
This motion, depending on the magnitude of its time scale, will be periodic or not . 

In the following we shall address the question of precession in accretion disks in systems 



similar to that considered by iRomero et all (120001 ). In our formulation, however, we will be 



concerned to treat the secondary tidal torque on the disk as a perturbation, which means 
that the secondary never crosses the disk, being far away from any point in the disk. Granted 
this, we will take into account the precession contribution to the total angular momentum 
and will discuss the conditions under which we may take an average of the Euler equations, 
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in such a way as to keep algebraic solutions to the precessional motion. Deviations from the 
precessional motion will be treated by solving the equation for the time evolution of 8, the 
angle that measures the misalignment between the normal to the plane of the disk and the 
normal to the plane of the secondary orbit. The time evolution will be obtained assuming 
the time scale for this motion is large or comparable to the remaining time scales. We obtain 
results for the precession velocity for different disk models, taking into account deviations 
from geometrical thinness due to a relation between the angular velocity of the disk and the 
height scale. One of the main results of this work is to show that, if precession is to occur in 
the inner region of disks around a supermassive black hole, the disk is not keplerian, being 
very thick and luminosity deficient. 



2. The System and the Dynamics 

The system we consider is a binary one composed by a primary compact object (a black 
hole), of mass M p , a massive secondary object (another black hole), of mass M s , and an 
accretion disk around the primary one. The secondary is orbiting the primary plus disk in 
a plane inclined of an angle 8 relative to the normal of the plane of the disk. Working in a 
reference frame tied to the disk, a point on it is at a distance R from the center (the primary 
star) and has coordinates (x, y, z) or (r,<f),z). The line of nodes will be chosen coincident 
with the x-axis, and the position of the perturber (secondary) will be 

d = d (cosuti + cos8 sinu t j + sinS sinutk^j (1) 

where cu is the angular velocity of the secondary. 
The potential at the point (r, 0, z) is 

^-Z^-^AU (2) 
R \R-d\ d 3 ' 

G being the gravitational constant. The first term is the potential due to the primary, 
the second is the tidal potential due to the secondary, and the last takes into account the 
acceleration of the coordinate system. The tidal force per unit mass will be 

CM 2 

f r = -GM s A(r, z, t, 0, d, u, 5)- 3 / 2 C(r, d, t, u, 0, 8) - B(t, u, 0, 5), (3) 



- 6- 



U = GM a A(r, z, t, 0, d, oj, 5)- z ' 2 dQ(t, oj, 0, 5), (4) 

and 

GM 2 

f z = -GM s A(r, z, t, 0, d, oj, 5)~ 3/2 N(z, d, t, to, 5) — - E(t, u, 5) (5) 

Where A, B, C, E, N, and Q are given by 

A(r,z,t,(f),d,u,5) = (d 2 + r 2 + z 2 -2rdB(t,w,<f>,6) - 2 z d E(t,u,5)), (6) 

B(t, oj, 0, 5) = cosuj t cos(p + cos5 sinus t sin<f), (7) 

C(r,d,t,u,<j>,6) = r-dB(u,<j),6), (8) 

E(t, 5, oj, t) = sinS sinu t, (10) 

N(z,d,t,oj,5) = z - dE(t,5,ou,t), (11) 

Q(t, oj, 0, 5) = —sin<f> cosoj t + cos5 sinu t cos<p. (12) 

To obtain the the precession angular velocity we must take into account its contribution 
to the angular momentum. In Fig.l, we present a schematic representation of the system we 
are considering. 

where z is the coordinate normal to the plane of the disk and Q , the precession angular 
velocity, makes an angle 5 with z. Therefore, the angular momentum is no longer parallel to 
the normal to the plane of the disk, and we may write 

L — L\\ + L ± 



— hoTi + I 2 oj 2 , 



(13) 
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where 1\ and I 2 are respectively the moments of inertia relative to an axis normal to the 
plane of the disk and another that lies in the plane of the disk (z=0); cJ*i and cJ* 2 are angular 
velocities parallel to those axes. It is straightforward to see 

io\ = Qa + Q p cos5, (14) 

u 2 = fl p sin5, (15) 
fid being the angular velocity in the disk. 

So, applying these considerations to an element of matter ( a ring) in the disk, we obtain 



dT — Q p A dL 



= d p A cJl d h + % A dfadh (16) 



or 



\dT±\ = Q p fid sinS dl\ + fl p sin5 cos5 (d I\ — d I 2 ) , (17) 



with d r± being the differential element of torque perpendicular to the angular momen- 
tum, d I\ and d I 2 , the moments of inertia of the ring. 

To obtain the precession velocity for a ring, it suffices to solve only for an element of 
torque perpendicular to z, say t x . Clearly, 

dh = pr 2 dV 
dl 2 = p (r 2 sin 2 (f) + z 2 ) d V. 
dr x = p {—z (fy cos(f> + f r sin<f>) + r f z sin<f>) d V. (18) 



Finally, we may write 
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sinS cosS (r 2 cos 2 (j) — z 2 ) Vl p 2 + r 2 Vl p Vtd sin5 + z (/^ coscj) + f r sintfi) — r f z sincj) = 0. (19) 
A rapid inspection in this equation tell us that: 

a- there is a constraint on(r, 0, z) that should be satisfied in order to have a solution for 
the precession velocity. If this constraint is satisfied we may assume the angular precession 
velocity approximately constant in time, as long as the time associated to its derivative is 
much longer than the precession period. If these conditions are not met, we should solve the 
differential Euler equations; 

b- if these constraints are satisfied, the precessing disk will admit double solutions. 
These solutions will be functions of (r, 0, z, t, to). 



It should be remarked that this treatment differs, e.g., from that of iKatz et all (119821 ). 



Papaloizou and Terqueml (119951 ) in the sense that 



- it takes into account the precessing angular velocity contribution to the angular mo- 
mentum of the disk, - it obtains constraints, conditions to the validity of such a treatment, 
- and the solution is double valued. 



3. The Constraints On The Disk 

From equation ( 19 ) , of the previous section, one can easily see that the constraints 
on the precessing disk comes from the inequality 



r 4 sin 2 8 Q d 2 + 4 sin5 cos5{r 2 cos 2 (f) - z 2 ) t x > 0. (20) 
Setting z = £, where £, the height scale of the disk , is related to p, the density, through 

8 7T \ Ik l 6 

with Qk being the keplerian velocity, M is the accretion rate, S a function that takes 
into account boundary conditions on the angular momentum and on the torque, and is given 
by 



S = 1 - fa (- 

\ r 



ri\V 2 
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/3 being the angular momentum in units of the keplerian one, (3q being this ratio at the 
inner radius (0< (5< 1), the inequality may be worked out to read, approximately, 



(3 2 r sinS M p d 2 p-Ar K 2/3 sinS (M p + M s ) (cosS ( 1 + 3 cosS T - + 3 sinS (^j j - 1 j p 1/3 



-12Kcos 3 5- 7 > 0, (22) 
a 

and 



(3 2 r d 2 M p sinS p-A cos 2 5 r 2 (M p + M s ) K 1/3 p 2/3 + A cos 2 5 (M p + M s ) K > 0, (23) 
where 

if = 7.510 29 1 pT - g. (24) 
a 

The solution to the first inequality is, approximately, 



r JM V + M S \ K , s 

p > 12cosfcoM (^_j_ ( 25) 



Concerning the second inequality, if 



(16) 2 / r V 2r 2 fM p + M s \ 2 , 

» cos 2 Scot 2 5 — ^— — - < 1, (26) 



27 VW V M p 

it will hold everywhere in the precessing disk. Therefore, using the definition of p, we 
may rewrite 

,.> 1MmM (M*)fA' (27) 



and 
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Since r > £, the solution is given by the last expression. Otherwise, if 

(16) 2 / r \ 4 9r ojM r) + M s \ 2 , x 

-^T\— J co^cot 2 ^ p n/f a ) >1, (29) 

we obtain 



27 VW V M p 



c JM„ + M S \ K K 



or 



(3 2 d? ~ r ~ r 



which may be rewritten 

— (^) G) 3 * ^ * ivi— ©' <*> 

It should be remarked that we may change inequality ( 32 ) into a condition on the 
masses, which reads 

1.25/? 2 rfi 5 2 < M P (M P + M S ) cos 5 cot 5 < 0,11 f.f 15 , (33) 

Zd M p 

where Zd is the size of the disk in units of Schwarzschild radius and the distance between 
the primary and the secondary is expressed in units of 10 15 cm. 



4. The Average Effect of the Density 

Our main interest are systems with r << d, a situation in which we may treat the 
tidal torque as a perturbation. Under this condition we may expand A(r, z, t, 0, d, u, 5) 3 ^ 2 , 
given by (eq.(6), to second order in r/d and z/d, and since the constraints obtained in the 
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previous section are satisfied, we may take the average over <fi and t on eq.(21) . After some 
straightforward but tedious algebraic manipulations, we obtain 

sin 5 cos 5 ^ — z 2 ^j Q p 2 + r 2 sin 5 f3 Q k Q P + u s 2 cos 5 ^z 2 — r =0, (34) 
where u s is the keplerian velocity of the secondary. 

We now specialize this equation for the average on z and write it in terms of the density 
p, i.e., 



p 2/3 r 2 ^6 cos5 tt 2 + 12 (3 tt k tt p - 3 u s 2 cotb = 2 K 2/3 (2 cos5 a/ - 3 cu 2 cotd) . (35) 
To see the behavior of fl p in terms of the density, let us first define the parameter p as 

1 27 \l3dj \ M r ) K 1 

and to have a better insight into the precession problem in accretion disks and to 
summarize the results we have obtained so far, we write the characteristic densities p c , p m in, P* 
and the actual density, p, in terms of p and po, the value of p for disks, i.e., 

K 

Pc = — , 



Pmin — 3.9 p I p c , 

a 



p* = 17.54 p 3 / 2 p c , 



P = \j) Pc=PoPc- (37) 

For p — > 0, there is no constraint and our treatment will be valid for any density. As p 
grows, but still < p < 1, there is a characteristic density in the problem and our treatment 
only holds for p > p m i n . For p > 1, there appear two more characteristic densities in the 
problem, p c and p*. In the region p c < p < p* our treatment doesn't apply as well. 
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For disks we are interested in, such as the a-standard , slim and advective models, 
1 < Po < oo. The equality at the left occurs when the disk gets thick either due to a huge 
accretion rate or advection. In that situation, p — > 1 

Now, we plot £l p = Qp(po), assuming M p >> M s , <5 4°, and d/r = 100. 

Figure 2 above shows that, in the upper branch, fl p decreases in the range < p < 0.1, 
and then increases in 0.1 < p < oo; in the other two branches, Q p grows with increasing den- 
sity. For p -> 12cos(fcofr^ M ^ p Ma ) ^gl^g, Q p ±cj s (| -^) 1/2 , respectively, in the upper 
and middle branches. For p — > | , f2 p — > ±oo in the lower and upper branches. In the middle 

branch, tl p - | (|f - 2) . Asymptotically, tl p - - (l ± (l + ff| cu s 2 coM V ^ . 

Angular velocities are in units of the keplerian one. 

Therefore, we will have, in the middle branch, 

'3 1 \ 1/2 „ 1 u v 2 cot5 d 

in the lower branch, 



U -\?7M £Si ' £ 4— P (38) 



^ ( ( cos5 o „ d\ . 



and, in the upper branch, 



/3 1 \ i/2 

£f! ' <0 °- (40) 

From equations (38) to (40), we can see a huge spread in the value of the precession 
velocity due to the effect of the density. 

Since for actual disksp = p p c and 0.01 < - < 1, the solutions are at the right of 
Po = 1. As a matter of fact, for the a standard model, the actual density is much greater 
than p c , the asymptotical behavior is a good approximation. Using this limit, we may see 
that differential precession is drastically reduced as compared to Katz 1973 result. 

It should be remarked that the constraints we have are much more restrictive than those 
obtained directly from equation (35), which are , approximately, 



p > 0.356p 3/4 p c - 



(41) 



13 



Now, we ask ourselves if, within this formalism, it is possible to recover iKatz et all (119821 ) 
result for the precession velocity, i.e., 



n 



K 



3 av 

"4 a 



cosS. 



(42) 



To answer this question we insert the above expression for the precession velocity into 
equation (35), to obtain 



Po I — I ^- I cos 3 S — 93 cos5 — 3 cotS — | = - ( ) cos 3 S — 6 cotS. 
1 8 VfiJ H r J 4 Vn fc ' 



(43) 



Since << it implies po ~ ^ , or 



2\/2 



1.5 



(44) 



Therefore, we must conclude that, under our formalism, IKatz et all (119821 ) results only 
applies to disks with density profile quite different from that expected under the keplerian a 
thin disk approximation. As a matter of fact, the results are in complete disagreement with 
this approximation. 



5. The Importance of deviations from the mean precessional motion 



Differently from IKatz et all (119821 ). in ou formalism, deviations from the mean preces- 
sional motion can't be treated as perturbations, i.e., nodding motions. As will be shown, the 
amplitude of these motions, though small for the precessional time scale, is quite large, as 
far as the appropriate time scale is considered. As a matter of fact, these deviations may not 
be periodic motions, but motions in which misalignment increases with time. The prevailing 
scenario will depend on the time scales involved. 

The importance of the misalignment motion in our formulation can be seen in the time 
evolution of 5. From equation (15), 



5 = Q p sin5. (45) 

Now, suppose that for some unespecified reason, Q p is given by the extreme right of the 
middle branch. Then, the time evolution will be given by the solution of 



-14- 



*=lm c ° s6 -r' (46) 

which is 

sin5 = y PQkr I , (47) 

where B is an integration constant. From Figure 1, we see that 5 is the angle between 
z and z', respectively the normal to the plane of the disk and the normal to the plane of the 
secondary orbit. 002 lies in the plane zz', and is normal to z. Therefore, if 5o is the angle 
between z and z'at t = 0, we must have, at any time, 



cosS = cosS cosQ 2 , (48) 

where G2 is the angle in a plane that contains z and is normal to 002 . Assuming that 
for t — > 0, 6 2 — * 0, we obtain 

B = (1 + sin 2 5 ± 2 sinSo) , (49) 

COS z 0r 



and 



_ .. B^ 2 expbt . 



cos5 V ^ ex P 2 6 1 + 1 / ' 



where 



1 a;, 2 d 



2 r 

We see that, as time goes on, the misalignment increases. 

Now, if we take our solution from the lower branch, assuming the unperturbed disk to 
be quasi keplerian, 



5 = — 2 p flk tanS, 



(51) 
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and 



sinS = Bexp2[3Q k t. (52) 

Proceeding in the same way as we did before, imposing the same boundary conditions, 
we obtain B 2 = sin 2 S and 



cos 2 6 2 = — -^r (l - sin 2 5 exp 2(3 Q k t). (53) 

COS 2 0q v ' 



It is easily seen that this solution holds for t < t , where 

*° = ~WFT in \sin5 \, (54) 
and for (5 ~ 1, in a very fast way, 62 — > 0.5 7r. 

Finally, let us analyse the situation (5 very small in the middle branch. This implies 

( d \ 1/2 

{277^5 ) ' (55) 

and 

5 = A{sin5) l/ \ (56) 
where A is defined in the previous equation. Making the substitution 

cos5o cos&2 = x, (57) 

yields 

A. (58) 



(l-x 2 f 4 ' 

Since we are only interested in the behavior of the solution, we make the approximation 

<( i - x2 )> 1,, (i4^j" A ' (59) 
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where <> stands for an average. Then, 

1 /Cexp-2 6t-l 

cosB 2 = 



cosSo \C exp — 2 bt + 1 J ' 
1 + cosSq 



C 



cos5n 



< {1-x 2 ) >V4- 
Again, the solution only holds for t < t , where 



(60) 



t * , n (i±£££*y (61) 



2 6 \ 1 — cos^c 

To obtain the time evolution of 5 we have assumed that the time scales associated to it 
are much longer than the remaining time scales in the problem, so any function of S survives 
to the averaging process. In other words, even after the averaging process, for any/ = f(5), 
5 is considered the instantaneous value. However, if the time scales are comparable, equation 
(45) should be interpreted as 

5 =<Q p sin5 >, (62) 

the solution being 



5 =< Q v sin5 > t + Sn 



(63) 



and 



cos02 = cos(< Q p sinS > t) — tg5o sin(< Q p sin5 > t). (64) 

Clearly, we have not considered viscosity in our formalism. Its inclusion certainly in- 
validates the use of eq.(45) to study the temporal evolution of the inclination angle 5. We 
should use, instead, a second order differential equation. 
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6. How Keplerian Thin Disks Precess? 

Most of the work on the literature assume an keplerian disk precessing under the perturb- 
ing influence of a third body, and to answer that question we shall take a naive approach to 
the angular momentum transport in the disk. If j in , j and j out are, respectively, the inwards 
angular momentum transport ( per unit time), the rate at which the angular momentum 
flows into the central compact object and the outwards flux of angular momentum, we have 



and we may write 



jin = JO + jout, (65) 

';,„, S M VTTJh: (66) 
jo = A) M y/GMn (67) 



and 



j out = 2i\r tr r4n (68) 
where T r( p is the stress tensor given by 

r r4> = -2r}a rt<t> , (69) 
ay t< j, being the rate of strain tensor, given by 

rd(3V k 

a ^ = _ 2^-' (70) 
where 14 is the keplerian velocity. 

Then, the angular momentum conservation may be written as 

(3MVGMr = (3 M ^GMr 1 + 27ir£T r4> . (71) 



In the above equation, we have assumed null boundary condition for the torque at the 
inner radius r — r\. This yields for the energy dissipation, neglecting (3 dependence on r, 
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3 GMM 2 ( A p (n\ l l 2 



^Z^r^-jW )• (72) 

and for the density 

3 Mp ( p /ny/A 

Now, we write the energy conservation, as matter in the disk moves from a point r 2 to 
a point r, as 



AE P + AE k + A£ int + Q" = 0, (74) 

where AE P is the variation of potential energy, AEk the variation of kinetic energy, 
AE int is the variation of the internal energy and Q~ is the total amount of energy that 
leaves the system from r 2 to r. Writing for the radial velocity 

V r = -a(^j V k , (75) 
we obtain, assuming r 2 » r, and that the heat produced by shear leaves the system, 



• GM MV k 2 ( o /A 2 2 {£V\ 3 o2 .GM/ Po fn\ 1/2 \ n ,„s 
- M — + ^r-[^ + {r) +a {r) ) + 2^ M — { 1 -j{t) J^ ' (76) 

where, again, we have neglected P dependence on r . Collecting like terms, we get 

4^-2A/>(^) ,/, + g)' + a»(g 4 -2 = 0, (77) 

which tells us that the keplerian approximation only holds for r = 1 and a 2 (^) 4 << I. 
Under the assumption of thin disk and a « 1, for r — > oo, /? — >~ 0.7. 

Now, assuming /3 constant, and using inequality (29), i.e., 

;32 <_ cos , OT( ^^_j(_), (78 ) 
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leads to 



d\ 2 64 cosScotS (M p + M s . 

Z) <~ 7^772 77774 ( Mp h (79) 



r) ~ 3^ 2- (^) 2 -a 2 (^) 4 
and, since d » r , this implies 



, / n\ 2 

,2 



or, approximately, since a « 1 , £ ~ r, which means the breakdown of the keplerian and 
thin disk approximations. Despite not having treated rigorously the angular momentum 
transport, this result is quite general, as it may be seen by considering the r-component of 
the force equation, 



1/2 ^!_^_I^_™ (80) 
or r p or r z 



which, after using the expression for p when the disk gets thick, is integrated to give 



3 a 2 - 



+ 5[- 

r 



6(1- /3 2 ) =0, 



^1) 



yielding 



(;) 2 =4(i-3 2 )- (82) 

To analyse the situation p < 1, the sign of inequality (78) is inverted, i.e., 

W -3v / 3 2-(^) 2 -^(^)H M p ^ 

and since d/r >> 1, a trivial solution is £/r << I. Therefore, we are led to the 
conclusion that p > 1 implies /? << 1, and p<l,/3«l. 

it should be reminded that the inner region of accretion disks around supermassive 
black holes gets thicker either due to the radiation pressure dominance or to the ineffective 
radiative cooling, a situation in which advection becomes the dominant cooling mechanism, 
and f3 < 1. 
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From these results, and those of the previous section, we must conclude that advective 
disks will precess with a velocity close to that given by p ~ p c , i.e., 

Qp * -db " { p2 + 1/6 ^ 2 coM (§7 " 2 ) ) V2 1 ' (84) 

and the standard-a keplerian thin disk will precess, roughly, with the asymptotic value 
of the precessing velocity, i.e., 



7. The p- profile of the accretion disk or precession and disk model 

The rate of precession depends on the kind of disk model we are assuming. Different 
models give different relations between p and po, as a consequence different precession rates. 
Besides, we have seen (3 is related to the scale height of the disk through 



P 2 = l- 



which means 



/3 2 = i-Rr 2/3 . (86) 

Inserting this relation into the definition of p, equation (62), to rewrite it as 

27 V A7p J (?o 2/3 -l) 

For the a-standard model, assuming black-body emission and gas pressure dominance 



Ma 



3/8 



Po w 4.53 x 10 11 ( ) » 1, (88) 



and 
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256 2r 2r /M„ + M s \ /n 4 

Mg and Mi are, respectively, the mass of the primary in units of 10 9 solar masses and 
the accretion rate in units of 1 solar mass per year. 

Therefore, under these assumptions, the disk is keplerian and our formalism holds every- 
where. The precession velocity will be given by the asymptotic limit of the middle branch, 
with (3 = 1, 

^p-J^f- cot6 - ( 9 °) 

To finalize, we can easily conclude that the a-thin keplerian disk is not affected by the 
constraints on the density. This, certainly, is not the case if the disk is thick. Let us now 
examine a disk radiative pressure dominated. Let us also assume cooling due to radiation. 
In that case, 



p 

ri 



P r = ^iF r , 
cm H 

P= l -p£l k H 2 . (91) 
This implies, using equation (72 ), approximately, 



r 8ircm H r \ \r 




the solution being 



r 2 



L((l+4a 2 ) 1/2 -l), (93) 



or 
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where 



Po 1/3 = t 2a i/2 , , 04) 

((l + 4a 2 ) 1/2 -lj 



0.34 Mi , x 

a= —w,- (95 > 



and z, now, is the radial distance in units of 3 R s . 

In figure 3, below, we plot p = Po(a). 

Using equation (87 ), and imposing p = p*, we obtain 



p 4 / 3 - 2p 2 / 3 + 1 = 4 1 4.45 cos 2 5cot 2 5 (Me±M±\ QV /3 , (96) 
which, using equation (94 ) above, changes to 

((1 + 4a 2 ) 1/2 - I)' -8a 2 ((l + 4a 2 ) 1/2 - l)' + 16a 4 = 

6631.2 cos 2 5 cotH ( Mn±Ml) V-^V a 4 . (97) 



v M p J \d 15/ 
Expressing in terms of the precession period, T p , and z in terms of a, yields 

fzV _ 1.95 x IP' 5 M[ / p 2 / 3 - 1 

UisJ " cos 2 5T p 2 a 3 (M p + M s ) 2 M p 3 \ Po 2 /3 

Inserting into equation (97 ) and defining 



(98) 



3.13 x 1(T 3 cos 2 5M? 
T 2 M 9 £ 



!/ = ^YT7- 5 L ' (") 



we get 



2/ = a(4a 2 -((l + 4a 2 ) 1/2 -l) 2 ). 



(100) 
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Below, a plot of y — y(a). 

A glance at this figure shows that if we are looking for a solution with p ps 1, we should 
choose around a ps 25.74. This gives 



p ss 1.06, 
Mi 

z d « 0.0132 — i, 

Mg 



di 5 ~ 5.07 



\ 1/4 

M x 4 (M p + M s f cos 2 5 cot 2 5 



M p 6 



M 3 \ V2 

T p = 2.81 x 10~ 3 cosS ( | , (101) 
where Zd is the size of the disk in units of 3 x R s . 

It is worth remarking that an acceptable solution gives z > 1, or ^ > 75.76 if the 
primary is a Schwarzschild black hole, or ^ > 25.25 for an extremely rotating Kerr black 
hole. Assuming these constraints do apply, we may write 

\ 4 

Z 



p ~ 0.15 ^— J ,l<z<z d . (102) 

However, if they don't apply, p > p* everywhere. Then, use of equations (89 ) and( 98 
) will lead to 

3 2 

p = 2.99 10- 3 cos 2 5 -4— 2 -— . (103) 

T - Ms Ua 2 -((l + 4a 2 f 2 -l) 2 J 

Since p given above has to be smaller than p given by equation (102 ), we obtain a lower 
bound on T p , given by 

T p 2 = 0.02 cos 2 5 ^4 "7 (104) 

p M 9 2 / 



^4 a 2 - ((l+4a 2 ) 1/2 -l)^ 
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Below, we plot T p 2 for 5 « 4, z d = 100 and M 9 = 6.0. 

Finally, if the disk is radiation pressure dominated and cooled by advection, equation 
(92 ) changes to 

a(^j + (1 -2.2a) (J-^J -a = 0, (105) 

the solution being 

- = — ({I - 4.4 a + 8.84 a 2 ) 1/2 - (l-2.2a)Y (106) 
r 2 a V y / 



or 



po i/3 = (l-4.4a + 8.84a 2 ) 1/2 + (l-2.2a) ^ {m) 



2a 

and equation ( 100 ) changes to 



y = a 3 ^(1 - 4.4 a + 8.84 a 2 ) 1/2 + (1 - 2.2 a)) 2 - 4a 2 ^)/ 

((1 -4.4a + 8.84a 2 ) 1/2 + (1- 2.2a)) 2 , (108) 

where y and a have been previously defined. 
Below, we plot p = po(a) and y = y(a). 

Proceeding as before, we look for a solution with p « 1. This will force us to look 
around a 0.4525, yielding 



Po ~ 1-01, 
/3 « 0.1, 
y w 5 x 10~ 3 , 
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T„ « 0.626 coso" 



. x V2 

M 



5 



2\ V4 

l s) 

M fl 6 



di 5 « 10.38 cos 2 8 cot 2 5 — — g-^- . (109) 



Reasoning; as we did before, & > 1.33 if the primary is a Schwarzschild black hole and 



Mg 

Mi 



> 0.22 for an extremely rotating Kerr blak hole. The p profile in the disk is still given 
by equation (102). Again, if these constraints on the density don't apply, p will be 



p = 4.7 10- 3 cos 2 5 -4—^ v — i: : -, (110) 

l p Mg 



3 ((1 -4.4a + 8.84a 2 ) 1/2 + 1 -2.2a) 

^(1 - 4.4a + 8.84a 2 ) 1/2 + 1 - 2.2a) 2 - 4a< 

and the lower bound on the precession period 

( (1 - 4.4 a + 8.84 a 2 ) 1/2 + 1-2.2 a) 
T p 2 = 0.031cos 2 o" ^ ^ -. (Ill) 

!/r (7(1- 4.4a + 8.84a 2 ) 1/2 + 1 -2.2a) -4a 



Figure 5 gives T p = T p (a) for 5 « 4, z d = 100 and Mg = 6.0. 



8. Conclusions 



We have treated the problem of precession in accretion disks taking into account its 
contribution to the total angular momentum. We have looked for conditions under which 
the problem may be treated by solving the algebraic Euler equations, i.e., constant angular 
velocities. We have found that the problem is characterized by the parameter p, given by 



256 cos 2 5 cot 2 5 [M„ + M s \ 2 fr\± 
p= ' ' 



27 p 4 V M v 
and the densities 



r\ 4 
d) ' 



Pmin = l^COsScotS 



+ M S \ K 



\ Mp J(3 2 d?' 
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K 



h\2 cos 6 5cot 6 5 



M„ + M,W r \ 6 K 



M p J \l3dj r 3 



For p < 1, the only constraint is p > p m i n , and the solutions, as a function of the density, 
lie in three branches: the upper one, with prograde precession velocities, at the upper left 
of p = |p c ; the lower branch, with retrograde velocities, at the lower right of p = |p c ; and 
the middle branch, between the upper and the lower, in which the solution is retrograde for 
p < ^Pc, and prograde for p > ^Pc- For p > 1, part of the p-space, p c < p < p*, is not 
allowed. Using simple energy and angular momentum transport arguments, together with 
constraints obtained in that formulation, it is shown that, for p > 1, a — > 1 and £ — > r, 
leading to the breakdown of the keplerian and thin disk approximations. In that situation 
the disk is luminosity defficient. As a matter of fact, a very important contribution of this 
paper is to show the incompatibility between keplerian disk and thick disk. The procedure 
we have adopted is suitable for treating deviations from the mean precessional motion, the 
misalignment motion, as well. For different limiting expressions for the precession velocity, 
we have obtained the time evolution for 5, the angle that measures the misalignement of 
the plane of the disk and the orbital plane of the secondary star, the misalignment angle. 
It has been shown that the time evolution of 5 depends on the magnitude of the time scale 
associated to this motion. Periodic motions only occur if this time scale is comparable to 
the remaining ones. Otherwise, there will be a tendency to force the plane of the disk to 
become normal to the secondary orbital plane. It is shown that, if the disk is keplerian, full 
misalignment, i.e. 5 = 0.5 n, is reached in a time comparable to the precession period, which 
means that the solution for the system is no longer given algebraically. We should make 
resort to the solutions of the differential Euler equations, a formidable task far beyond our 
goal in this paper. Finally, for different disk models we have shown that the constraints on 
the density will hold as long as ^ >> 1. In that situation, we have found the properties of 
the precessing region as a function of Mi and M 9 . An expression for the separation distance 
between the primary and the secondary is also found. Otherwise, if ^ < 1, lower bounds 
on the precession period and on the separation distance are found. These bounds depend on 
jj^. Application of this formalism to different disk models, may be summarized as follows: 
1- for the keplerian thin a standard disk model, the precession period (in years) will be 

T p > 0.08 z^Mp, 
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For accretion disks in binary systems, we may take the size of the disk as the truncation 
radius, the point where the Roche equipotentials first intersect. In our systems we are not 
allowed to do so but, in any case, we found in the literature that values for z ~ 100 — 1000 
are not unusual. Besides, expecting precessiuon periods of about hundred years, the system 
will reach full misalignment in a time of the order of the precession period. The disk is not 
affected by constraints on the density. 

2- For disk in which radiation dominates pressure and cooling, if -jg 1 > 75.76 for a 

Schwarzschild black hole and H 1 > 12.62 for a Kerr black hole, density in the disk may come 
close to the critical density p*. In that situation, t pa r, and the size of the disk can not 
exceed 



z d = 0.0132 ^. 

Mg 

Besides, if the system is not to be affected by gravitational radiation, should be 
much larger than the values given above. Precession period will be 



T„ = 2.81 x 10" 3 cosS 



Mf 



1/2 



Mt 



5 



A reasonable value for Zd would imply a very large ratio We are, therefore, led to 
the suspicion that the assumption I pa r, or po pa 1, is very strong for this model. Abandon 
of this assumption leads to a lower bound on the precession period, given by 



T p 2 = 0.02 cos 2 5 ^4 

M 9 2 



^4 a 2 - ((l+4a 2 ) 1/2 -l)^ 



3-Finally, for a radiation pressure and advective cooling dominated disk, I pa r, and its 

size 



Mi 

Again, if the system is not to be affected by gravitational radiation, otherwise it will 
live for a short period of time, ^ >> I. If ^ is not too large, £ « r, and the constraints 



on the density will not be effective. Again, there will be a lower bound on the precession 
period given by 



T 2 = 0.031cos 



((1 - 4.4 a + 8.84 a 2 ) 1 ' 2 + 1-2.2 a) 



Mi 



( 




-4a 2 



) 



To finalize, we would stress the main contribution of this work as pointing out the 
possibility of considering the disk structure, when studying disk precession, by means of an 
alternative procedure. 



Caproni, A., and Abraham, Z. 2002, submitted 
Katz, J.I. 1973, Nature, 246, 87 
Katz, J.I. 1980, ApJ, 236, L127 

Katz, J.I., Anderson, S.F., Margon, B., and Grandi, S. 1982, ApJ, 260, 780 
Katz, J.I. 1997, ApJ, 478, 527 

Larwood, J.D., Nelson. R.P., Papaloizou, J.C.B., and Terquem, C. 1996, MNRAS, 282, 597 

Larwood, J.D., and Papaloizou, J.C.B. 1997, MNRAS, 285, 288 

Larwood, J.D.1997, MNRAS, 290, 490 

Livio, M. 1999, Phys. Rep., 311, 225 

Lubow, S.H.1991, ApJ, 381, 259 

Lubow, S.H. 1992, ApJ, 401, 317 

Lubow, S.H. 1992, ApJ, 398, 525 

Lubow, S.H., and Pringle, J.E. 1993, ApJ, 409, 360 

Ogilvie, G.I. 2002, MNRAS, 330, 937 

Ogilvie, G.I. 2002, MNRAS, 331, 1053 



REFERENCES 



-29- 

Papaloizou, J.C.B., and Pringle, J.E. 1977, MNRAS, 181, 441 
Papaloizou,J.C.B., and Terquem, C. 1995, MNRAS, 274, 987 
Priedhorsky, W.C., and Holt, S.S. 1987, Ap&SS, 45 , 291 
Pringle, J.E. 1996, MNRAS, 281, 357 
Pringle, J.E. 1997, MNRAS, 292, 136 
Roberts, W.J. 1974, ApJ, 187, 575 

Romero, G.E., Chajet, L., Abraham, Z., and Fan, J.H. 2000, A&A, 360, 57 
Vishniac, E.T., and Diamond, P. 1992, ApJ, 347, 435 



This preprint was prepared with the AAS IATgX macros v5.2. 



-30- 




Fig. 1. — Binary System Geometry. P stands for primary and S for secondary; zand z' are, 
respectively, normal to the plane of the disk and to the secondary orbital plane; 5 is the 
inclination angle between them; d is the secondary distance to the primary. The line of 
nodes is coincident with x 
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Fig. 2. — Solution for the precession velocity, displaying the three branches. Under the 
assumption of hydrostatic equilibrium, only the middle and the lower branches are allowed 
( Po > 1 )■ We have assumed M p » M s , 5^4, and d/r =100 
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Q. 




Fig. 3. — po = (^) as a function of a = ^ j± for a Schwarzschild Black hole (assumed in 
the plot), and a = ^ ^ for a maximally rotating Kerr Black Hole. 
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Fig. 4. — y, defined by eq.(lOO) and eq.(108), as a function of a = ^ assuming the 
primary as a Schwarzshild Black Hole 
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Fig. 5. — The precession period squared ( year 2 ) as a function of a = ^ |^-, assuming the 
primary as a Schwarzshild Black Hole 



